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ABSTRACT 

The winter-time transport of geostrophically computed 
angular momentum at 500 mb for the period 1 December 1955 
to 28 February 1956 is obtained for Sault Ste. Marie, 
Michigan. This transport is subjected to a cospectral 
analysis to compute the contribution of harmonic waves with 
periods 1 - 30 days to the total angular momentum transport. 
Computations are made for two other stations, one east, and 
one south of Sault Ste. Marie for comparison of transport as 
related to general circulation processes, 

The author is deeply indebted to Professor F. L. Martin, 
Department of Meteorology, U.S. Naval Postgraduate School, 
for his encouragement and invaluable aid in the preparation 


of this manuscript. 
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Symbol Definition 


AQ increment equal to one centimeter of length along a 
parallel of latitude 


St) mean cross covariance of two continuous functions 
of time 
oy lag applied to a continuous function 


U(t) time series of the wind component along a parallel 
of latitude 


w(t) time series of the wind component along a meridian 
i averaging period 

gu (nh) Fourier transform of u(t) 

Gar (n) complex conjugate of the Fourier transform of ar (t) 
nN harmonic number 


é. cospectral estimate 
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Definition 
frequency 
chi~square value at the 5% level 


mean cospectral value between peak and trough of 
the cospectrum 


degrees of freedom 


Nyquist frequency 
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1. Intreduction 

Widger [10] in his study of the flow of absolute 
angular momentum in the atmosphere stated: 

Although even a complete knowledge of the angular 

momentum trausfer in the atmosphere cannot by itself 

furnish a solution of the problem of the general 
circulation, any proposed scheme for this circulation 
should include a means for securing the angular 
momentum flows which are observed. For this reason, 

it is essential that all possible observational knowl- 

edge concerning these processes be obtained. 

Many authors in recent years, Widger,‘Starr, and 
Mintz [10, 7, 3] for example, have provided studies in 
connection with the relationship of the absolute angular 
momentum transport to the general circulation, or the global 
babanee of aengular.momentum Inegenmeral thesewestudies 
provide quantitative findings of the spatial distribution of 
the total momentum transport and establish the fact that 
there is a mean poleward transport. 

This paper, in hopes of COMIN ne to the knowledge 
of angular momentum transport, has as its purpose a cospect- 
tral analysis of the angular momentum transport Ber waves of 
periods 1 - 30 days. Due to the limited time available for 
the study, it was confined to the 500-mb level at a station 
in the mid-latitudes. The primary station employed was 
sault Ste. Marie, Michigan. 

Several other cospectra are included for other 
stations for the purpose of comparison. However, the 
spectrum for these comparative stations are for periods 


a 


< - 30 days. 
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2. Mathematical Formulations 

Widger [10] has shown that the rate at which the ab- 
solute angular momentum of the zonal wind is increasing in 
a portion of the atmosphere bounded by the earth's surface, 
the top of the atmosphere and extending from the pole to a 


given latitude $ may be expressed as 


ae ; f 
at JPMdIV ES PAU dS +f prdot{[ruds 


where dV is any element of volume within the polar cap just 
described,dS an element of the surface of the volume,do& 
ime projection of Migeteooe of the earth's surface area on 
the meridional plane passing through a given point, Pp is 
the pressure, r the radius of the earth, and Ya the east- 
ward frictional stress "acts only om the air ay vhe cama. 
muriace. 

In the first term on the right V, is the meridional 
component w of the velocity into the polar cap. Since the 
vertical velocities at the earth's surface and top of the 
atmosphere may both be considered zero, this quantity as a 
whole represents the rate at which angular momentum is 
brought into the region by air motiom across the boundary. 
The second and third terms on the right represent, respect-— 
ively, the rate at which angular momentum is changed due to 
the torque exerted by pressure differentials across mountain 
ranges and the westward torque,rT, , exerted by frictional 
drag at the earth's surface. lIowever, since this study 


will consider only the transport through a unit area of the 





latituce wall at the 500-mb surface, the terms involving 
pressure and friction torques may be neglected. No attempt 
is made herein to determine the left side of equation (1) 
from the computed transport at 500 mb. 

The absolute angular momentum of a gram of air at 


latitude ¢ may be written 
; A. y 
y\ = er cbs ow 4 cos* J 2) 


Furthermore, since V,=w and the element of surface may be 
taken as an elemeut of the latitude wall, dSsdgdz2, taking 
the first term on the, right in equation (1), the transport 
of absolute angular nomentum may be written as the sum of 


GMOs terms, 


roos ff fae da dz (3) 


equal to the rate of advection of "relative" angular 


momentum, and the du*momentum defined as 
ae. ricos*g ff par AA az (4) 


which is equal to the rate of advection of angular momentum 
due to the earth's rotation. Since ft-momentum does not 
hend»itself to a cospectral analysiseand will not be cem- 
puted in this study, the chief emphasis is placed upon the 
average angular momentum transport through a unit area and 
per unit time given by equation (43). Therefore, after 
applying the geostrophic approximation for U andW in 
equation (3) and dividing by N, one obtains the average 


angular momentum transport 
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where the hydrostatic equation has also been used in (3). 
At first glance equation (5) appears to be a summation | 
of N terms around the latitude circle at a distance AX 
apart, and thus a representation of equation (3). However, 
a slightly different interpretation was placed upon equa- 
tion (5), in that AX represents a grid distance of 300 
nautical miles centered on the station for which the 
computation is being made, and 4% and Az, are the contour- 


height increments defining the mean geostrophic ‘wind com- 


ponents 
U.2-7 ae and ty er : (6) 


in the neighborhood of the station. In addition, AY was 
taken equal-to AX in equation (6). Thus the interpretation 
of equation (5) actually employed here was to regard it as 
a time-average, over N equal time intervals, of values of 
relative angular momentum transport centered at the station. 
The vertical increment for convenience was taken as AP equal 
to one millibar and AA on the right side of equation (5). 
was scaled down to one centimeter instead of the equivalent 
length ‘of 400° -nautreal mi besewsedMior the sritd pe lfictmics 
of angular momentum transport in this study will then be 
Cun? nab . 

Equation (5) shows that the geostrophic relative mom- 


entum transport depends upon the cross covariance between 
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U3 and WV. It therefore represents the average geostrophic 
relative angular momentum transported across a section of 
the latitude wall with srid length AAs Jars and vertical 
thickness Ap= / mb. 

The mean cross covariance S(t) of two stationary, 
Gaussian time-series W(t) and a(t), at lagt, may be ex- 


pressed by the convolution integral 
res a ber p 
S(r) - 7H é ff u(t) ar (tet) dt (7) 
O 


Frequently in meteorological applications, S(t) is 


"evenized" by forming: the quantity 


Sin+ Ser. (8) 
2 


which is independent of the sign of lagV¥. Henceforth, it 
will be understood that the mean cross covariance has been 
@venized. 

The cious (or simply cospectrun) - of the 
two series u(t) and ar(t)is defined as that part of the mean 
cros¢ covariance associated with a unit wave band centered 
at wave numbern in both series. In symbols oe may then 


be defined as 


-_ 
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where [' is the averaging period and 9y(n) and gar(n) are 
the Fourier transforms of u(t) and w(t) respectively, de- 


fined as in 


u ~A0Qt 
9,2 43e fTuaye*@* dt 9 we 2Z2 20) 
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The circumflex indicates the complex conjugate. 

If one now employs the definition of cross covari- 
ance Sev) for an evenized function, as contained in equation 
(7), the Faltung theorem leads to the Fourier integral 


expression 
oe T 2 RY 
200 Sy 


so that S(w) and Gan are spoken of as Fourier-transform 
pairs. Equation (11) may be inverted to solve for Cy in 
terms of Sir). Also, the resulting inversion cam be spec-— 
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1alized to the case of a finite, but long, averaging period 
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A proof of equation (12) is given, from first principles, 
imeAppendix is 

When discrete data are employed at time intervals 
of At , using a maximum lag of mat, the maximum period, 


or fundamental, that can be resolved is 


Twn ~ Gn Et GR 
In equation (12), since TV appears as a summand, Yt will be 
written as 


v=hat gal Ve | (14) 


with #& taking on each of the discrete values shown in 
(14). Since the averaging period of equations (12) is 


le mat, equation (12) leads to 
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fenem-el 
In the case of Ge and Coy , the resulting coefiiciemys are 
divided by two, since only one-half bandwidth is available 
at m=O and nz7m. Combining equations (5), (7), and (8), 


the cross covariance for k lags is 


My 2 Fn cee D2 AZ yp + AZ, cag, AZ, ; | AA AP (16) 
Oofkem 
since the averaging period of equation (7) is now =(N-k) At, 
The factor two in the denominator of equation (16) results 
from evenizing. 

The cospectral estimates offered by equation (15) 
are termed "raw estimates". Blackman and Tucey [1] have 
iseecued several smoothing procedures, the simplest of 
which involves application of a movingeweighted average def- 


red iy 
Ci Se Coe Pe, a + Se Ge Ca? ) 


In the case of N=M there is no higher value to include in 


the average and the weighted average is modified so that 


~~. 


Comat alownci! Haque (18) 





4. General Procedures 

The initial step required establishing what data 
were to be used in making computations. Recalling that the 
computed periods are determined by the relation Zm At/n - 
where M™ is the maximum number of lags at discrete time 
imtervalerAct , it its clear that @ Half-day Pac interval re- 
quires 40 lags to obtain a spectrum for periods 1 - 40 days. 
Cn the other hand, a one-day sample and lag interval re- 
quires only 15 lags to obtain a spectrum for periods 2 — 40 
days. However, Blackman and Tukey [1 | state that the best 
results are obtained when the longest lag is about 5 - 10 
per cent of the length of the record. This would require 
that the retord, if possible, consist of 300 -— 600 observ- 
ations when using the 500-mb charts produced at 0300 GMT 
and 1500 GMT daily, or half this number when using only one 
Series of «charts. in order tovebtamn this number vem 
observations, it would require a minimum of 5 months data 
and the spectrum would be influenced by seasonal regimes. 
In order to reduce the data-processing load and inter-annual 
influences as much as possible, it was decided touse the 
eeerratwens from the three winter months of December 1955, 
January 1956, and February 1956. Thus the longest lag was 
only 17 per cent of the length of the fecord for Waits 
particular season. However, no claim is made that the 
cospectra obtained are necessarily representative of those 
for other winter seasons. 

Data were taken from two sources, depending on the 


station being considered. since heights of both the 0300 
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GMT and 1500 GMT 500-mb level were required at Sault Ste. 
Marie, Michigan, the primary station studied in this paper, 
facsimile charts transmitted by the U. S. Weather Bureau 
wervemmralizéd. In the case of the other stations used for 
comparison, the smallest period computed was only two days 
and the data used were the 1500 GMT charts published in 
Daily Series Synoptic Weather Maps [a]. 

From the above-mentioned charts, the geopotential 
height differences were obtained with a 5.0 degree (lati- 
tude) grid. Assuming that normal synoptic smoothing was 
applied in preparation of the charts, a grid of this dimen- 
Sion should not filter out significant short waves that 
have not already been removed by synoptic smoothing. 

In five cases, where the 0300 GMT facsimile charts 
were missing, geopotential height differences were computed 
from the actual observed wind listed in Northern Hemisphere 
Data Tabulations [9]. The observed winds were reduced to 
geostrophic values by subjectively compensating for the 
degree of cyclonic or anticyclonic curvature associated with 
the flow. 

The smoothed cospectra of geostrophic relative 
angular momentum transport were obtained by using these 
data in equations (15), (16), (17), and (18), and are graph- 
ically displayed in figs. 3 through 6. ~All computations 
except for the final smoothing, utilizing equations (17), (18), 


were accomplished on the CDC - 1604 digital computer. 
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4, Limitations and Sources of Error 

The computational results of this study have been based 
on discrete data using equal intervals At. In one series 
el Computations, indicated by the solid line dn fies 9, che 
data interval was reduced to 12 hours in order to determine 
the effects of sampling fluctuations and of atmospheric 
oscillations. In all other cospectral analyses, the data 
interval was 24 hours. Throughout the study the maximum lag 
was 15 days. 

The period associated with the smoothed cospectral 
estimate is 2mat/m with nN=I,2,-°°°: ,ma. On the other hand 
the central frequency f, of the frequency band is 

> nae ne Lamy 
aad the highest frequency resolvable occurs for naam j4 this 
frequency is called the Nyquist frequency i ; wae 


cospectral estimate corresponds to the frequency band 


(n~ ha) Jamar to (nt yf 2 mAS 

It should be not@d@ thigtcsaG @espdifficult to resolve 
frequencies contributing to the cospectrum at the longer 
periods. For example when h=/ , the frequency band spans 
the range !/gn cycle / day to /fo cycle / day. 

In order to assess the reality of sharp peaks or 
troughs, Panofsky’s criterion ior r/» distribution of the 
power spectrum has been followed. Strictly speaking a much 
more rigorous method of assessing significance would have 


been achieved by testing the cospectral correlation 
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Heme - 1s the cospectral estimate for the yn th harmonic , 
and En 8 and Pp, arare the power spectral estimates for t& 
and ~u- at the same harmonic. However, the power spectra 
computations were not available in time for the completion of 
this aspect of the work. Hence, the more qualitative ap- 
proach indicated in fig. 1 was attempted with some misgiv- 
ings. It must be emphasized that when the term "957% 
confidence limit" is seGip it is only by analogy with corres- 
ponding phraseology to be found in Panofsky and Brier [4 , 
Pp. 145 | »- The actual confidence may be at a somewhat lower 
level. 

Applying the XA distribution to the cospectrum, a 
peak or trough covld not have happened randomly (at the 95% 
level of belief) provided the peak or trough in the cospectral 


estimate lies outside the range AP, defined as 


AP:(P)(L/ ) (19) 


where FP is the mean cospectral intensity between the peak 
and trough, e y is the chi-square’ vailtie at tite .05) Bevel 
for “ “deerems of freedom (V =2N/m-VYo = }2) . Accord- 
ing to Panofsky and Brier [a, ion 55 | the limiting value of 
Raw? = 
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Fig. 1. Illustration of the X/v test of 


Significance at the 95% level of belief. 
A second source of error which may considerably bias 
the accuracy of cospectral estimates is a phenomenon called 
"aliasing". This phenomenon removes energy from higher 
frequencies than those sampled, and may convert it to lower 


frequencies. Fig. e illustrates how aliasing may occur. 





Tllustration of the effects of sampling 
variations in giving an aliased spectrun. 
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suppose a periodic function, such as that illustrated 
im es. 2, Exists for both Ut) ana ar(é) , and consider the 
period of the oscillation shown to be 2Mt. If lags are 
taken at intervals of 34€ , one is investigating the 
cospectral estimate at a period of EAs. ait ake sampling 
times corresponded to points A, B, and C, one obtains a 
rather large "aliased"cospectral estimate for the period GAT. 
Howevem, af the samplime pointe had ween By Fyrand G, the 
aimesed cospeetral contributioneatepeniod GAT would be 
zerg. Note that the suesnanigwen could be a real atmos- 
pheric oscillation or.meteorological "noise" of an obs- 
ervational character. The primary atmospheric oscillations 
of periods shorter than one day are the semi-diurnal (12-hr) 
wave and a 10.5-hr wave, originally described by Pekeris [5] : 
However, the amplitude of these waves is quite small. 

Blackman and Tukey [1] have shown that the frequencies 
which are the principal aliases of frequency a are Fi A 
Ut t fF, \ GFA e » eters, wiere Fin= M2At is the Nyquist freaq- 
uency. However, the effect of aliasing is most pronounced 
at frequencies fe near the Nyquist frequency. Griffith et al 
[2 suggest that one way of eliminating aliasing from the 
cospectrum is to delete the cospectral estimates for the 
last 20 to 40 per cent of the higher frequencies. However, 
an experiment was performed, using data for Sault Ste. Marie, 
to assess the extent of aliasing. 

Fig. 4 shows three cospectra for Sault Ste. Marie 
based on: 


(a) the 1500 GMT reports at 24-hr intervals (curve A), 
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(b) the 0300 GMT reports at 24-hr intervals (curve B), 
and 


(c) both the 0300 GMT and 1500 GMT reports at 12-hr 
imFervals (enmipue 6% 


It may be noted from fig. 3 that all cospectral est- 
imates at periods greater than 2.74 days are in fairly close 
agreement. However, curves 4 and B deviate rather markedly 
at their Nyquist frequency of % cycle / day. 

In the ¢ase of curve C there’ siowle bemiaaov le: om™no 
aliasing at periods greater than two days according to 
Griffith's criterion, since 50 per ‘cent of the spectrum 
lies between the periods one and two days. On the other 
ned , “atewall other statrons i wrllebe assumed, hence! omm , 
that all cospectral estimates at periods 2.73 days or larger 


are relatively unaffected by aliasing. 
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pee interpretation of *Cospectra 

Figs. 3 through 6 present the six cospectra computed 

for this study. In these figures the ordinate is the co- 
pectrum and the abscissa, plotted for convenience on a 
logarithmic scale, is the period. 

In assessing the reality of the various peaks observed 
in fig. 3, the procedure embodied in equation (19) will be 
applied only to curve ©, as the other two curves do not 
Gaeratve srgnifieantlywexcept ata period "oie, 75 days, which 
will be discussed separately. In order to compute a mean 
cospectral estimate fer the peak at a period of 10 days, it 
was assumed that the cosvectral estimate at 40 days is the 
apex of a trough. iis neenuet of thes 95%"and 5% fiducial 
limits established from the mean cospectral values on each 
side of the peak and trough, indicate that the peaks and 
trough associated with periods 10 and 6 days, respectively, 
are Significant. The peaks and troughs have only a 5% 
probability of randomly exceeding these limits; therefore 
both are ceomsadered to beerea ls 

In the case of curve A, the peak appears not to be 
Significantly different when ompared to the fidticial limats 
computed for curve C. However, limits computed from mean 
values of curve A, although not shown here in order to avoid 
confusion with those of curve 0, show the value at the peak 
on curve B is considerably greater than the limit to the 
left side and slightly greater than the limit to the right 
Side. Although the peaks at 4.75 days are not markedly sig- 


nificant from the limit on the right side, there is 
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justification for stating these peaks seem to be real when 
one considers both the left and right side limits. Using 
the significance test discussed above, the reality of any 
manor pealgs.in fag. vat periods of 4S daystor leeS\is 
doubtful. 

The curve in fig. 4 is the average cospectrum for 
Caribou during the period 1 December 1955 to 28 February 
1956. A one-day interval between data points was used to 
compute this cospectrum. An inspection of the fiducial 
limits shown in fig. 4 indicates that the peak at 6 days is 
Significantly different from the limits on both the left and 
right side of the peak, and it is therefore considered to be 
real. As to the small peak at approximately 3.3 days, the 
confidence interval shown to the left is quite representative 
for this peak. It can therefore be concluded that this peak 
(at 3.3 days) is probably not real, as it does not lie beyond 
the confidence range. However, it can not be inferred from 
this that the peak is a result of sampling variations, since 
the test does not with certainty disprove the reality of the 
peak. 

The curve in fig. 5 is also for Caribou, but was conm- 
puted from data taken during the winter previous to that 
used in fig. 4. The limits in this case show the troughwand 
peak at 7.5 ana 5 days, respecvIvely , * oMbemrertrrrre 07 
at 4.75 days is questionable, since it is not significantly 
different from the limit to the right. 

One of the more striking features in each of the 


cospectra is the large transport associated with the longest 
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period and the fact this transport has the same sign as the 
total transport. 

inmeeenésal itis dificult to -associate a particiuea: 
period of the cospectrum with a specific space wave, since 
the period is a function of both the wavelength and speed 
om thie Space wave. It, thereby, seems possible for a trans— 
poetwat a specific period to have contributions from both 
long and short space waves. This would be the case if the 
ratio of their wavelengths is equal to the ratio of their 
speeds. However, assuming that the Rossby wave equation 
gives a valid relationship between the wave speed and wave- 
length for any given mean zonal wind, it may be expected 
that long space waves move eastward slowly and short ones 
mome rapidly. In view of this assumption, the transport 
associated with a period of thirty days may be attributed 
to the contributions from long space waves. 

Furthermore, large values of transport by the longest 
period were to be expected, for an examination of the 
monthly mean 500-mb charts of December 1955 through 1956 
showed Sault Ste. Marie and Caribou were on the western 
side of a mean long wave trough, hence the net transport 
would most likely be southward (negative). However, for the 
second cospectrum at Caribou (fig. 5), the period December 
1954 through February 1955 was selected by ascertaining 
through monthly mean 500-mb charts that Caribou was east of 
the trough and a northward (positive) transport could be 
expected. In the case of Athens, Georgia, the main reason 


for selecting this station to provide a comparative 
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cospe¢trum is the fact it is almost due south of Sault Ste. 
Marie. One could expect these two stations to have a net 
transport of the same sign, but examination of the monthly 
mean 50Q-mb charts show the mean long wave trough to have a 
southwest to northeast tilt of sufficient magnitude to place 
the trough west of Athens. As a result the net transport for 
Athens was northward, while that of Sault Ste. Marie was 
southward. The indicated convergence of momentum transports 
Should then act in such a way as to increase the mean zonal 
westerlies. Note also that the only pronounced peak in the 
cospectrum at Athensioccurswvat a period of 0 das, and that 
the cospectral values diminish with decreasing period. ‘This 
is probably due to the fact that short period waves do not 
penetrate appreciably into this latitude belt. 
intberprebatwon of wagnificamt tramsports associated 
with shorter periods requires reference to the earlier 
statement, that short waves are likely to move rapidly. In 
view of this concept the momentum transport evident near 
periods of four days is attributed to waves of approximate 
space-wave numbers 6 to 10, which Saltzmann and Fleisher 
6] have called "cyclone waves". However, the momentum 
transport peaks at periods near 10 days are presumably due 


to long waves of space-wave numbers 1 to 5. 
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6. Vemelusions 

It may be concluded, tentatively, from the cospectra 
presented in figs. 4 through 6 and the preceding discussion 
CIA : 

(a) Maximum transport of geostrovhic relative angular 
momentum is accomplished by disturbances with the longer 
periods. 

(bd) The peaks with periods near 10 days are related 
to longs snace waves. 

(c) The direction of the net transport is highly de- 
pendent on the position of the mean long wave trough 
relative to the observing station. 

(d) Cyclone waves of space-wave numbers 6 to 10 may 
contribute to the maximum momentum transports associated 
with periods near 4 days. 

(e) The maximum transports at periods of 40 days may 


be attributed to the circulation index cycle. 
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Fig. 4 The cospectrum of the mean relative angular momentum transport at the 500-mb level 


over Sault Ste. Marie, Michigan for the period 1 December 1955 through 28 February 


1956. 
days. 


The auxiliary graph has the ordinate enlarged 100 times for periods 1 to 2 
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Fig. 4 The cospectrum of the mean relative angular momentum 
transport at the 500-mb level over Caribou, Maine for 
the period 1 December 1955 through 28 February 1956. 





























































































































Fig. 5 The cospectrum of the mean relative angular momentum 
transport at the 500-mb level over Caribou, Maine 


for the period 1 December 1954 through 28 February 1955. 
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Fig. 6 The cospectrum of the mean relative angular momentum 
transport at the 500-mb level over Athens, Georgia 
for the period 1 December 1955 through 28 February 1956. 
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APPENDIX I 
Derivation of the Cospectral Function 
Mie cress @ovarrance at lagt, of two functions cf 


time, is given by 
if 
ze | ; 
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The first cosine term integrates to zero for all cases 
m,n 4 0. The second cosine term integrates to zero for 


all cases m An. When m =n, integration gives 
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The evenized cross covariance equals one-half the sum 


of the cross covariance at ¥ and-T , 


cos (¥ - On) caso" An. Anzh2--- 
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Defining the unsmoothed cospectrum as 
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the Fourier transform of GG) is obtained as 


Sof Beet dh 


2 
Multiplying both sides by cos-~Z2t sent dv 


and integrating from 


Otol gives 


ef St) Co —— Seer dt = - en J Cth) Gag Sone) dv 


4Tr nt 
For all n # 0,CoS —-—— integrates to zero, giving 


, | eae er 
Cin= 2 al S(v) cos ERE dr 


where |’ is the integration limit and 7 is the fundamental 


werTod. 
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